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ABSTRACT.—Fifty-seven individual plasma samples from desert tortoises (Xerobates agassizi) representing 10 
separate populations were analyzed by polyacrylamide gel electrophoresis using alkaline buffers. An albumin-like 
protein was found to be polymorphie for two eleetromorphs in northern populations inhabiting the Mohave Desert 
Provinee, while Sonoran Desert populations to the south were monomorphic. The genetic divergence demonstrated in 
this survey is similar to earlier studies and provides evidence for the Colorado River as a potential barrier to gene flow 
among tortoise populations. These data suggest that tortoise plasma, examined by various electrophoretic methods, 
may provide a nondestructive means of determining the broad regional origin of desert tortoises. 


Desert tortoises (Xerobates agassizi) pres- 
ently inhabit two regions of southwestern 
Utah, separated east and west by the Beaver 
Dam Mountains. The population dynamics 
of the Beaver Dam Slope tortoises (west of 
the mountains) have been severely impacted 
by both human and animal activities during 
the past several decades. Despite both fed- 
eral and state protective regulations, tortoise 
numbers in the slope region are probably at an 
all-time low (Mike Coffeen, UDWR, personal 
communication). Two stable populations near 
St. George, east of the Beaver Dam Moun- 
tains, face new human development projects 
that threaten the future of these previously 
isolated populations. The eastern populations 
are found in Paradise and City Creek canyons, 
and relocation of some tortoises from these 
populations is presently under consideration 
by the Utah Division of Wildlife Resources. 
However, several relocation issues remain 
unresolved, including methods of collection, 
conditioning, transport, sex ratios, and speci- 
men numbers involved. One problem to be 


addressed is the question of genetie compati- 
bility of Utah's tortoise populations, especially 
if any of the tortoises east of the Beaver Dam 
Mountains are translocated to the western slope. 

Few data are available comparing the 
physiology or morphology of Utah’s separate 
tortoise populations. Rainboth et al. (1989) 
electrophoretically analyzed whole blood 
homogenates of 146 desert tortoises from two 
separate localities in California for allozvme 
expression at 23 loci. The two populations 
were quite similar, as each locality contained 
unique elements only when allozyme combi- 
nations were used and there was a consider- 
able degree of overlap in allozyme frequency. 
Lamb et al. (1989) included five tortoises from 
Paradise Canyon in an analysis of phyloge- 
ographic patterns in mitochondrial DNA 
(mtDNA) of the gopher tortoise complex. 
Their data indicated that these individuals fit 
within an eastern Mohave “elone” represent- 
ing the northern Arizona and eastern Nevada 
region. Also, Jennings (1985) included five 
tortoises from the Beaver Dam Slope region 
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TABLE 1. Genotypes and variability estimates for the polymorphic GP-1 locus resolved in 10 localities of Xerobates 








agassizt. 
Number of GP-] 

Locality N genotypes A* He 
12. Pima Co" AZ 3 3AA 1.0 0.0 
2. Tucson, AZ 6 6AA 1.0 0.0 
3. Pinal Co., AZ 4 4AA 1.0 0.0 
4. Maricopa Co., AZ 4 4AA 1.0 0.0 
5. Kingman, AZ 3 ZAA (BB 3 0.0 
6. Beaver Dam Slope, AZ 12 9AA:2AB:1BB 1S 0.083 
7. Paradise Canyon, UT 12 3AA:6AB:3BB NaS: 0.208 
8. Lincoln Co., NV 9 LAA: LAB Jat 0.250 
9. Riverside Co., CA 4 LAA:3AB 1B} O38 

10. San Bernardino Co., CA i SAA:4AB 15 0.214 





* — mean number of alleles/locus. 
** — mean heterozygosity-direct count. 


in biogeographic investigations comparing 
blood and tissue enzymes using horizontal 
starch gel electrophoresis. His study included 
no specimens from east of the Beaver Dam 
Mountains in Utah but showed that a north- 
to-south variation was evident and that the 
Arizona Beaver Dam Slope specimens fit 
within the northern (Mohavean) group. 

Mitochondrial DNA and isozyme studies 
are costly and may involve traumatic (to the 
tortoise) and labor-intensive biopsies of inter- 
nal tissue or the sacrificing of specimens for 
necropsy. The purpose of this study was to 
determine whether general proteins in tor- 
toise plasma could be used to detect geo- 
graphical differences in tortoise populations. 
Our primary interest was to compare plasma 
proteins among Utah’s aforementioned tor- 
toise populations using alkaline polyacry- 
lamide gel electrophoresis. However, as elec- 
trophoretic profiles of plasma or serum from 
Xerobates agassizi have not been previously 
reported, we examined plasma profiles from 
tortoises from Arizona, California, and Ne- 
vada as well. This report follows the taxo- 
nomie grouping of gopher tortoises by Bram- 
ble (1982) and the nomenclature revision of 
Bour and Dubois (1984), applying Xerobates 
as the genus for the desert tortoise. 


MATERIALS AND METHODS 


Study Area and Sampling 


A total of 68 desert tortoises were collected 
from 11 localities in Arizona (AZ), California 
(CA), Nevada (NV), and Utah (UT), repre- 
senting most of the range of this species in the 
United States. Plasma samples were collected 


in heparinized tubes (3 ml) by venapuncture 
(jugular vein or antebrachial sinus) from 15 
tortoises from the Beaver Dam Slope, AZ, and 
Paradise Canyon, UT. Additional plasma sam- 
ples were donated by colleagues involved 
in desert tortoise projects in other parts of 
AZ, and in CA and NV. Loealities and sample 
sizes for the populations sampled are listed in 
Table 1, and their geographie loeations are 
plotted in Fig. 1. 


Electrophoresis 


Plasma samples were lyophilized and stored 
at 4 C. Polyacrylamide gel electrophoresis 
(PAGE) was run in a Bio-Rad vertical-slab 
electrophoresis cell, Model 220. The gel was 
1.5 mm thick. Sample wells 10 mm long were 
formed with Canaleo stacking gel (2.5% acry- 
lamide). The 7% acrylamide separating gel 
was 80 mm in length. Samples were elec- 
trophoresed toward the anode (+) at 20 mA 
constant current and were stopped 10 mm 
from the bottom of the gel. The eleetrophore- 
sis buffer used was 0.025 M_ Tris/0.192 M 
glycine, pH 8.3; bromophenol blne was used 
as tracking dye. The gel was fixed in 10% 
acetic acid/40% isopropanol/50% water for 60 
min, stained with coomassie blue (0.05%) in 
10% acetie acid/10% isopropanol/80% water 
for 60 min, and destained in 10% acetic acid/ 
10% isopropanol/80% water, using three or 
four changes of solution over a 24-hr period. 

Only qualitative differences were exam- 
ined in this study, and only the faster migrat- 
ing proteins were used, since alkaline PAGE 
is not the preferred method for resolving 
differences in basic (globulin-like) proteins, 
due to their short migration distanees. The 
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Fig. 1. Distribution of collection localities of tortoise genotypes AA, BB, and AB listed in Table 1. Numbers indicate 


localities. 


electrophoretic mobilities of rattlesnake 
(Crotalus atrox) plasma albumin and bovine 
serum albumin were used as markers and 
compared with the albumin-like protein(s) of 
the desert tortoise samples. Both albumin 
markers migrated at faster rates than the 
tortoise albumin-like proteins analyzed in 
this study. 


Genetic Analysis 


Two protein loci, identified in order of de- 
creasing anodal mobility, were designated 
“general proteins (GP) ] and 2. Allelic data at 
both loci were recorded as individual geno- 
types for analysis with the BIOSYS-1 program 
of Swofford and Selander (1981). Measures of 
genetic variability computed for each popula- 
tion sample included average locus hetero- 


zygosity (H, direct count) and the mean num- 
ber of alleles per locus (A). The genetic dis- 
tance and similarity coefficients of Nei (1972, 
1978) and Rogers (1972) were calculated for all 
pairwise combinations of samples (corrected 
for small samples sizes as described by 
Levene [1949]), and all such matrices were 
clustered by the UPGMA algorithm of Sneath 
and Sokal (1973). Genotype ratios from the 
largest samples (localities 6 and 7, n = 12 for 
both) were tested for conformance to Hardy- 
Weinberg proportions by the X° goodness-of- 
fit option of BIOSYS, again corrected for small 
sample sizes (Levene 1949). 


RESULTS 


Sixty-eight individual samples were ana- 
lyzed from four states, including 11 localities 


4 J. L. GLENN ETAL. [Volume 50 


ihe 





GP- “28 
GP-1 => Ss . 





(+) 


AB 


AA AB A BB 
Fig. 2. Alkaline polyacrylamide gel illustrating GP-1 and GP-2 loci (enclosed in dotted rectangle), with allelic 
anion at GP-1 locus for plasma samples from three individuals (run in duplicate). Animals 1, 2, and 3 were 


consistently scored as genotypes AA, BB, and AB, respectively. 


BB 


TABLE 2. Matrix of genetic distance coefficients of Nei (1978, above diagonal) and Rogers (1972, below diagonal) for 
all pairwise combinations of Xerobates agassizi localities; locality numbers are as listed in Table | and D values are 
rounded off to two decimals. 





Locality ] 2 8 4 5 6 a 8 9 10 
] — 0.00 0.00 0.00 0.03 0.01 OA7 0.00 0.05 0.06 
2 0.00 — 0.00 0.00 0.03 0.01 0.17 0.00 0.05 0.06 
3 0.00 0.00 — 0.00 0.03 0.01 0.17 0.00 0.05 0.06 
4 0.00 0.00 0.00 a 0.03 0.0] 0.17 0.00 0.05 0.06 
5 GO1% 0.17 0.17 Oy a 0.00 0.00 0.00 0.00 0.00 
6 0.08 0.08 0.08 0.05 0.08 —_ 0.05 0.00 0.00 0.01 
7 O27 0.27 OF O20 0.10 0.19 _ 0.01 0.00 0.00 
§ Os 0.13 0.13 0.13 0.04 0.04 O15 — 0.00 0.00 
, 0.19 Oo19 LS, org 0.02 0.10 0.08 0.06 — 0.00 

10 0.18 O43 0.18 0.15 0.01 0.10 0.09 0.05 0.0] —_ 





(considered as separate populations). The pro- — df = 


tein profiles present in 7 of 9 plasma samples 
from the Desert Tortoise Natural Area (CA) 
and 2 samples from Utah were unique but 
were very likely artifacts due to their badly 
hemolyzed condition. These 11 samples were 
excluded from the data analysis. This reduced 
the number of sainples to 57 and the number 
of localities to 10 (Table 1). General protein-| 
was polymorphic for two electromorphs (des- 
ignated A and B) in several samples (Fig. 2). 
Table 1 summarizes the ratios of GP-1 geno- 
types across these 10 localities and the esti- 
mates of variability across both loci. Geno- 
type ratios at localities 6 and 7 conformed to 
Hardy-Weinberg expectations (X° 2.799) 
P = .095; X° = .503, P = .478, respectively, 


1 in both eases), suggesting that this is a 
simple Mendelian co-dominant system with 
two alleles segregating in some populations. 
Table 2 summarizes pairwise comparisons 
of two genetic distance coefficients (Rogers 
1972, Nei 1978) and shows that between- 
sample divergence was minimal. Neis D val- 
ues, for example, range from 0.00 to 0.17. 
Four of the Arizona samples (Maricopa Co., 
Pinal Co., Pima Co., and Tucson) are identieal 
(Neis D — 0.00). Five localities having the B 
allele at GP-1 formed a distinet separate clus- 
ter, albeit the total degree of divergence from 
the monomorphic populations was slight (D ~ 
0.05). Within this group, the Arizona and 
California populations were nearly identical, 
while the Paradise Canyon (UT) samples were 
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Rogers’ Distance 


Fig. 3. Dendrogram based on Rogers’ (1972) genetic distance values for 10 samples of Xerobates agassizi (see Table 1, 
Fig. 1). Clustering was by the UPGMA algorithm of Sneath and Sokal (1973), and the cophenetic correlation value was 


0.826. 


the most divergent. These relationships are 
also visually displayed in the UPGMA den- 
drogram presented in Figure 3, using the 
statistical analysis method of Rogers (1972). 
All other UPGMA dendrograms gave identi- 
cal or nearly identical topologies (data not 
shown). 


DISCUSSION 


Dessauer (1970) reported that while some 
blood proteins in reptiles are relatively 
conservative, others are quite polymorphic. 
Many species can be readily distinguished by 
the electrophoretic mobilities of blood pro- 
teins, and certain subspecies and populations 
can also be distinguished by differences in 
plasma albumin-like proteins (Dessauer and 
Fox 1958, Masat and Dessauer 1968). Masat 
and Dessauer (1968) also found that the 
albumin-like proteins of the Testudines have 
slower migration rates in alkaline buffers 


than do these same proteins in most other 


reptiles (and mammals). We also found that 
under alkaline conditions the desert tortoise 
albumin-like proteins migrate at slower rates 
than rattlesnake plasma albumin and bovine 
serum albumin. 

The results of this investigation suggest 
geographical differences in genetic variability 
of the albumin-like protein (GP-1) of desert 
tortoises. The northern (Mohavean) popula- 
tions were polymorphic, whereas the south- 


ern (Sonoran) populations were monomorphic 
at the GP-1 locus. An east-west Mohave dif- 
ference was observed due to the eastern isola- 
tion of the BB genotype (Table 1, Fig. 2) in 
populations from the eastern Mohave region 
of Utah and northwestern Arizona (Fig. 1). 
The B allele was not present in any of the 
central and southern Arizona samples. Those 
samples expressing the B allele may differ in 
the frequency of this allele, as suggested by 
the differentiation between Paradise Canyon 
and the Beaver Dam Slope, but the present 
sample sizes are too small for accurate deter- 
mination. Despite the smal] sample sizes, the 
heterozygosity estimates are more similar for 
the Paradise Canyon population and the three 
Mohavean populations (8, 9, and 10) than for 
the Paradise Canyon and Beaver Dam Slope 
populations (Table 1). If the Paradise Canyon 
tortoises differ from the Beaver Dam Slope 
population(s) in frequency of the B allele and 
are in fact more similar to California tortoises, 
this could reflect: (a) a divergence of allele 
frequencies between the slope and Paradise 
populations in allopatry, with allele frequen- 
cies at Paradise simply drifting to values 
similar to California populations (similarity 
by convergence); (b) transport of tortoises 
from California or Nevada to St. George and 
“dumping” into Paradise Canyon, but not at 
the slope (i.e., human-induced gene flow 
between California, Nevada, and Paradise); or 
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(c) transport and release of Arizona tortoises 
(most with AA genotypes at the GP-1 locus) on 
the Beaver Dam Slope but not Paradise 
Canyon, which would cause the A allele at the 
slope to increase in frequency at the expense 
of the B allele and “push” the slope population 
away from Paradise Canyon and California 
allele frequencies (see genotype differences in 
Table 1). These are not mutually exclusive 
hypotheses since a certain amount of “mixing 
may have occurred at all localities of Califor- 
nia, Nevada, Beaver Dam Slope, and Paradise 
Canyon populations. 

The long history of human collection and 
translocation of desert tortoises between 
states constitutes a variable that could influ- 
ence the genetic structure of desert tortoise 
populations, especially when comparing al- 
lele frequencies between populations well 
known as captive release sites. Hundreds, 
perhaps thousands, of tortoises have been 
picked up along roadways and released in dif- 
ferent regions over the past several decades 
(Mike Coffeen and Eric Coombs, UDWR, 
personal communication). This logistical dis- 
placement of tortoises continues at present. 
Certain localities have been popular release 
sites, e.g., regions of southern California 
(Desert Tortoise Natural Area), Arizona 
(McDowell Mountain region near Phoenix), 
Nevada (near Las Vegas and recreation areas), 
and Utah (St. George). The most common 
avenues of translocation by motorists crossing 
the Mohave Desert are east to west and vice 
versa. Favorable habitats for tortoises exist 
from Washington County in southwestern 
Utah to southern California—a region heavily 
traveled over a major interstate highway for 
decades. In addition to the release of tortoises 
in this region by motorists and local citizenry, 
Utah’s Beaver Dam Slope population(s) have 
been the site of approximately 200 captive 
release tortoises, regulated by the Utah Divi- 
sion of Wildlife Resources since 1970 (Mike 
Coffeen, UDWR, personal communication). 

Although Utah's tortoise populations lo- 
cated east of the Beaver Dam Mountains are 
often regarded as “captive escapees in a 
nonindigenous setting, there are no scientific 
data confirming this view. The hypothesis 
seems to have originated from several 
sources, e.g., magazine and newspaper arti- 
cles, and the opinions of a few naturalists 
and herpetologists. Support for the intro- 
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duced status is based on three general obser- 
vations. One is the fact that some tortoises 
found in the St. George region obviously have 
been captive specimens, exhibiting rope or 
chain holes in their shells or having painted 
areas on their shells. Second, residents of 
Washington County often keep tortoises as 
pets, and some have escaped or were pur- 
posely released. Third, the Beaver Dam and 
Virgin Mountains presently form an east-west 
barrier between natural assemblages, and a 
few reptiles found west of the mountains are 
not present east of them (e.g., Crotalus scutu- 
latus, Phyllorhynchus decurtus, Dipsosaurus 
dorsalis). The possibility still remains that 
some of Utah's tortoises found east of the 
mountains may be derived from ancestral 
stock of naturally occurring tortoise popula- 
tions that have since mixed with captive re- 
leased specimens. Support for the natural 
population relies on the fact that the Mohave 
Desert Province extends into this region, and 
many other Sonoran life-zone animals found 
on the western slope are also found east of the 
mountains and are natural assemblages. Like 
the desert tortoise, some of the reptiles are 
life-zone specific and occur on both sides, for 
example, the banded gecko (Coleonyx varie- 
gatus), Gila monster (Heloderma suspectum), 
and sidewinder (Crotalus cerastes ). 

The geographical differences observed in 
this investigation are similar to those found in 
the allozyme survey by Jennings (1985), the 
mtDNA survey by Lamb et al. (1989), and the 
morphometric analvsis of tortoise remains by 
Weinstein and Berry (1987). Specifically, the 
present study supports the earlier molecular 
investigations of Lamb et al. (1989), which 
showed divergence between tortoise popula- 
tions north and west of the Colorado River and 
those to the south and east. Their report pro- 
vided good evidence that tortoise populations 
now isolated on opposite sides of the Colorado 
River have likely been separated from each 
other for several million years. The mtDNA 
lineages from central and southern Arizona 
formed a single haplotype that differed from 
the northern haplotypes in CA, NV, UT, and 
extreme northwestern AZ by a minimum of 17 
restriction site changes (see Fig. 2 in Lamb 
et al. [19S9]). This is one of the highest levels 
of intraspecific genetic divergence reported 
for any animal species and exceeds that 
reported for many interspecific comparisons. 


1990] 


The exception in our stnudy was the small 
sample from near Kingman, AZ (locality 5 in 
Fig. I), which genotypically grouped with the 
Mohavean populations north and west of the 
Colorado River. Consequently, the genotypic 
composition for population 5 must be inter- 
preted with caution. The single BB homo- 
zygote in a sample of three individuals would 
not be expected unless the B allele was segre- 
gating at a high frequency. These results may 
be due to several factors, e.g., the transloca- 
tion of this specimen by humans, a sampling 
error for a low-frequency allele, or a degrada- 
tional artifact in this sample. If future sam- 
pling verifies the presence of a_ high- 
frequency B allele at this locality, it could 
represent an ancestral polymorphism shared 
with Mohavean populations north of the Colo- 
rado River. This anomalous result under- 
scores the need for statistically adequate sam- 
ple sizes in all future genetic studies of the 
desert tortoise. Therefore, identifying the 
specific origin of any individual tortoise on 
the basis of nuclear gene markers may be 
difficult. Weinstein and Berry (1987) suggest 
using a combination of physiological and mor- 
phometric screening methods to designate re- 
gional types. They analyzed shell morphology 
of adult (> 180 mm) desert tortoises of both 
sexes by using morphometric data gathered 
by the Bureau of Land Management at 31 
different localities. These measurements were 
collected from tortoise remains by several 
persons over a 48-year period. These authors 
noted that live tortoises were not used in their 
analysis and that shell morphology of tortoise 
remains does incur some shrinkage over time 
following death. They recommend further 
studies comparing live tortoises. None of 
these authors suggested that the differences 
observed in their investigations justify sub- 
specific designation for any of the regional 
populations, and thus Xerobates agassizi re- 
mains a monotypic species. 

Alkaline PAGE was useful in examining the 
albumin-like proteins in tortoise plasma but 
does not resolve slight differences in the elec- 
trophoretic mobility among the majority of 
the plasma proteins. However, this method 
did detect the polymorphic nature of the albu- 
min-like (GP-1) protein, and this protein may 
be one “marker” that could be used for desig- 
nation of broad regional types. Since the re- 
sults of this investigation are, with the possi- 
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ble exception of locality 5 noted above, very 
similar to the findings of others in that the 
broad regional genotypes in desert tortoises 
are approximately concordant, the PAGE 
screening of the plasma protein marker may 
provide one inexpensive method of objec- 
tively determining the regional origin of 
tortoises. If allele frequency data are to be 
used, additional specimens from Paradise 
Canyon and Beaver Dam Slope populations 
are needed to determine the significance of 
the allele frequency differences between 
these two localities. Additional samples from 
throughout Nevada and California would also 
be required. Also, some variations observed 
in the slower-migrating proteins could be ex- 
amined with more high-resolution techniques 
(e.g., isoelectric focusing, two-dimensional 
electrophoresis) combined with bio-image an- 
alytical instrumentation that can quickly and 
accurately scan and record qualitative differ- 
ences in electrophoretic profiles. Morpho- 
metric data from live tortoises should be col- 
lected to compare Paradise Canyon and 
western slope populations in conjunction with 
future molecular analyses. In fact, external 
morphology (such as shell shape) could be a 
more functional criterion for the survivability 
of Paradise Canyon tortoises on the Beaver 
Dam Slope (see Weinstein and Berry [1987]), 
since genetic differences between these pop- 
ulations appear to be slight. 
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